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Abstract: In this report we deal with the problem of global output feedback stabilization of 
a class of n-dimensional nonlinear positive systems possessing a one-dimensional unknown, 
though measured, part. We first propose our main result, an output feedback control pro- 
cedure, taking advantage of measurements of the uncertain part, able to globally stabilize 
the system towards an adjustable equilibrium point in the interior of the positive orthant. 
Though quite general, this result is based on hypotheses that might be difficult to check in 
practice. Then in a second step, through a Theorem on a class of positive systems linking 
the existence of a strongly positive equilibrium to its global asymptotic stability, we propose 
other hypotheses for our main result to hold. These new hypotheses are more restrictive but 
much simpler to check. Some illustrative examples, highlighting both the potential complex 
open loop dynamics (multi-stability, limit cycle, chaos) of the considered systems and the 
interest of the control procedure, conclude this report. 
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Stabilisation globale d'une classe de systemes positifs 

partiellement connus 



Resume : Nous considerons ici un probleme de stabilisation globale par retour de sortie 
d'une classe de systemes non-lineaires positifs de dimension n possedant une partie de forme 
analytique inconnue, mais cependant mesuree. Nous proposons tout d'abord notre resultat 
principal, un controle par retour de sortie, qui stabilise globallement le systeme vers un 
equilibre reglable appartenant a l'interieur de l'orthant positif de l'espace d'etat. Bien 
qu'assez general, ce resultat est base sur des hypotheses qui peuvent se reveler difHciles a 
verifier sur un exemple concret. Dans un second temps, par l'intermediaire d'un theoreme sur 
une classe de systemes positifs liant existence et stabilite globale d'un equilibre fortement 
positif, nous proposons d'autres hypotheses permettant d'appliquer notre resultat. Ces 
nouvelles hypotheses, bien que plus restrictives, sont nettement plus faciles a verifier. Le 
rapport s'acheve par des exemples qui illustrent a la fois la potentielle complexity des 
comportements en boucle ouverte qui peuvent etre produits par les systemes consideres 
(multi-stabilite, cycle limite, chaos) et l'interet de la procedure de controle proposee. 

Mots-cles : Systemes positifs, systemes partiellement connus, stabilisation globale. 
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1 Introduction 

Positive systems of ordinary differential equations are systems producing trajectories, that, 
if initiated in the positive orthant of R™, remain in this orthant for all positive time. Such 
systems are being payed much attention since they abound in many applied areas such as life 
sciences, social sciences, chemical sciences, telecommunications, traffic flows etc... [ill 13. Ill] . 

In this paper we deal with the problem of global output feedback stabilization of a class 
of n-dimensional nonlinear positive systems possessing a 1-dimensional unknown, though 
measured, part. Our main result generalizes an earlier one obtained in lower dimensions 
and for simpler structures in the context of bioprocess control [l^ . 

This paper is organized as follows. We first introduce some notations and definitions 
related to positivity for vectors and dynamical systems. We then introduce the class of 
uncertain positive systems that is concerned with our control problem. We then propose an 
output feedback procedure, taking advantage of measurements of the uncertain part, able 
to globally stabilize the system towards an adjustable equilibrium point in the interior of 
the positive orthant. The hypotheses on which our main result is built might however be 
difficult to check in practice. In order to derive some hypotheses that are much simpler to 
check, we propose a theorem on a class of positive systems (the strongly positive concave co- 
operative systems) that links the existence of a positive equilibrium to its global asymptotic 
stability (GAS). Some examples found in the litterature illustrate the potential complex 
open loop dynamics that might be produced by the systems under study (multi-stability, 
cyclic behavior or even chaos). Simulations on periodic and chaotic behaviors stabilization 
illustrate the control efficiency. 



2 Preliminaries 
2.1 Notations 

We consider the following autonomous nonlinear dynamical systems in R n : 

i = /(*)• (1) 

For the sake of simplicity, we assume that the functions encountered throughout the 
paper are sufficiently smooth. 

In the sequel, we will use the following notations: 

• x(t, xq) denotes the forward orbit at time t of system initiated at xq. 
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• Df(x) denotes the Jacobian matrix of system (QJ at state x. 

• Cf denotes the Lie derivative operator along the vector field defined by system QJ. 

• £|xj=o denotes a state vector whose i— th component equals zero. 

Moreover, a matrix is said to be Hurwitz if all of its eigenvalues have negative real parts. 
2.2 Positivity for vectors and dynamical systems 

As we deal with positivity in R™ and some related concepts, we first clearly state what we 
refer to as positivity (resp. strong positivity) for vectors and dynamical systems. 

Definition 1 (Positivity) 

x G M. n is positive (resp. strongly positive) and denoted x > (resp. x 3> 0) iff all its 
components are non-negative (resp. positive). 

We denote H.™ and int(M") the sets of positive vectors and strongly positive vectors 
belonging to K", respectively. Same definitions of the relations > and ^> (i.e. component 
by component) will be used for matrices too. 

We now state the definition of positive systems of ordinary differential equations as well 
as the definition of a special class of positive systems, the strongly positive ones. 

Definition 2 (Positive Systems) 

System Q) is a positive system (resp. strongly positive) iff: 

Vi G [l..n], Vx^o > 0, Xi(x\ x . =0 ) = fi(x\ Xi = ) > 0, (resp. fi(x\ Xi=0 ) > 0). 

It is straightforward that Definition [2] implies the following: 

Proposition 1 

Consider a positive system £7J) (resp. strongly positive). Then: 

V(x ,t) G K" x M.+ (resp. V(x ,t) G x E+), x(t, x ) > (resp. » 0). 

Remark 1 Though we have equivalence between Definition^ and Proposition^in the case 
of positive systems (this is the original definition of positive systems by we only have an 

implication from Definition^ to Proposition^ for strongly positive systems. Indeed systems 
with fi(x\ Xi= o) > and a strongly connected Jacobian Matrix may produce strongly positive 
trajectories too Q/. We however keep Definition as the definition of strongly positive 
systems since this criterion is much easier to check on dynamical systems. 
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3 Main result 

3.1 The considered class of systems 

In the sequel we focus on the (global) asymptotic stabilization problem of the following class 
of nonlinear positive uncertain systems in R": 

x = uf(x) + ap(x), ,^ 
V = i>{x). 

• u > being the scalar input of system 

• / : R™ -> R n 

. ceM" 

• ip : R n —> R is an output of system (5J 

Function ip(.) features the uncertain/unknown part of the system. For control purposes, we 
suppose that, though analytically unknown, might be online measured. We moreover 
assume the following: 

Hypotheses 1 (HI) 

1. Vz e int(]R!Ji), ip(x) > and Vi, Ciip(x\ Xi=0 ) > 

2. /(.) is such that the system x = fix) is positive 

3. 3/3 TO G M + such that V/3 > (3 m , the system: 
x — [3f(x) + c. 

is a strongly positive system and possesses an equilibrium x*^ which is GAS on 
R™ . 

These hypotheses will be commented after the statement of our main result. For now, 
just notice that Hypotheses Hl-1 and Hl-2 together with u > ensure the positivity of 
the considered class of systems (2J. 

3.2 Main result 

Theorem 1 

Consider a system 0) and let Hypotheses HI hold. Then for all 7 > f3 m , the nonlinear 
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control law: 

u = jy = jt/j(x). (3) 

globally stabilizes system \2ty on int(R™ ) towards the strongly positive equilibrium x* such 
that: 

fix*) = — c. 

1 



To prove Theorem we need the following Lemma on Lyapunov functions for positive 
systems possessing a strongly positive GAS equilibrium. 



Lemma 1 

A positive system possesses a strongly positive equilibrium point x* GAS on int(R" ) iff 
there exists a smooth real valued function V(x) such that: 



Vx » 0,x ^ x*,V{x) > 0, 
V(x*) = 0, 

Vx » 0,1 ^ x*,V(x) < 0, 
\J{x, V(x) < a} = int{Rl). 

a>0 



Proof (Lemma IT]) 

Consider a positive system and suppose it possesses a strongly positive equilibrium x* 
GAS on int(R J |). 

Now consider the following increasing (in the 3> sense) change of variables: 



In 



(x*) 



In 



/ 



that maps the strongly positive orthant on R™ and x* to 0. Then we have: 

x* 



(4) 



Since x* is a GAS equilibrium on int(R") for so is on R™ for system Then, from 
Kurzweil's converse Lyapunov theorem on GAStability (see Theorem 7 in there exists 

a smooth radially unbounded Lyapunov function W{z) for system iQJ. Coming back to the 
original variables, consider the function: 



V(x) = W(ln(x) - ln(x*)). 
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Since W(z) is a Lyapunov function for system ( |4 |. we have: 

Vx e int(R'l), a; ^ a;*, V(x) = W(z) > 0, 
y(a;*) = W{0) = 0, 

V.t e int(R^), a; ^ a;*, F(x) = W(z) < 0. 

Moreover, since W(z) is smooth so does V(x) and since W-^z) is radially unbounded on R", 
we have: 

\J{z, W(z)<[3} = M. n . 

/3>0 

Thus: 

|J {x, V(x) < «} = int(R^). 

which concludes the first part of the proof. 

The reverse implication is easily obtained through classical Lyapunov theory. □ 

Proof (Theorem IT]l 

Consider the positive system {3$ under control law ©■ We get: 

x = tp(x)(jf(x) + c) = g(x). (5) 

Since 7 > (3 m we know from Hl-3 that the system: 

x = -ff(x) + c = h(x). (6) 

is strongly positive and possesses a (strongly positive) GAS equilibrium x*. Then, from 
Lemma m for all 7 > /3 m , system © admits a real valued function V 7 (x) verifying Lemma 
properties. 

Now we consider system (0. From Hl-1, it is clear that we have: 

Vx » 0,x ^ x*, V 7 (x) > 0, 
V-y(x*) = 0. 

Moreover: 

CgVy(x) = 1p{x)C h V 1 (x). 

so that since Hl-1 holds: 

Vx>0,x^x;, C g Yj(x) <0. 

and the following still holds: 

\J{x, Vy(x) < a} = intQR"). 

a>0 

So that we can apply Lemma to the controlled system Q. Then the equilibrium x* is 
GAS for the controlled system © on int(R+). □ 
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Remark 2 It is important to notice that the proposed control law does not, to be applied, 
require any a priori analytical or quantitative knowledge of the function ip(x). The only 
two requirements are (i) its positivity (a qualitative property), and (ii) the possibility to 
measure this quantity on-line. This is particularly important for instance in biological process 
control, when some parts of the model are only qualitatively known (usually in terms of sign) 
while some other parts are precisely (analytically) known. The special case of bioprocesses 
presented by [ity does illustrate this point and shows the interest of such control procedures 
for real life bioprocesses management. We come back to this point later in the examples 
section. 

3.3 Comments on Hypotheses HI 

As we have previously noted, Hypotheses Hl-1 and Hl-2 guarantees the positivity of the 
systems 

First part of Hypothesis Hl-1 defines a qualitative property of the unknown part of the 
system, which, though one-dimensional, possibly acts on all state variables dynamics via 
the vector c. Notice that this is a really loose hypothesis since the only required property 
is ip(x) positivity as x is strongly positive. No specific analytical form is assumed on 
The second part of Hl-1 is required to ensure system (|2jl's positivity as the input u equals 
0. Hl-2 is required for the system's positivity as the input u is positive. 

Hypothesis Hl-3 is a stronger one. It reads that if the uncertain part tp(.) were positive 
and constant, then there must exist an input infimum /3 m above which system (|2J would be 
strongly positive and would possess a (strongly positive) GAS equilibrium. This property 
might recall "minimum phase" conditions, frequenlty encountered when dealing with the 
stabilization of nonlinear and/or uncertain systems. Minimum phase conditions ensures the 
GAStability of an equilibrium as the system's outputs are constant. Minimum phase based 
control results usually stabilises the output (s) in order to conclude on the whole system, using 
results on cascaded nonlinear systems' stability. However uncertainties for these systems 
are in general not on the outputs' dynamics, but rather in the intrinsic "stable" part. A 
"minimum phase" approach can thus not be used in our case since the output ip(x) and its 
dynamics is, to a large extent, unknown. 

From the analytical expression of the system that is to be stabilized, both Hypotheses 
Hl-1 and Hl-2 are easy to check. This is however not the case for Hl-3. The existence of a 
unique and strongly positive equilibrium of: 

x = /3f(x) + c 

might indeed be easily proven, but the demonstration of its global stability is in general a 
much harder task. Nevertheless, for some special classes of positive systems, one can link 
the existence of a strongly positive equilibrium to its global asymptotic stability. This is 
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for instance the case for positive linear systems, see e.g. jlll . |j|. It is clear that such a 
property is of great interest: it allows to prove in one step both the existence of a strongly 
positive equilibrium and its GAStability i.e. Hypothesis Hl-3. Some works can be found 
in the litterature regarding this problem for nonlinear systems: see flftl These results 
can indeed be used to check Hl-3. It is to be noticed that these works are all related to 
cooperative systems (see Definition 0] in the sequel). 

In the next section we also base our approach on cooperative systems. We propose a 
class of positive systems for which the existence of a strongly positive equilibrium implies 
its GAStability. We then derive a set of "easy to check" hypotheses that are sufficient for 
our control procedure © to be applied. 



4 Sufficient conditions to verify HI 
4.1 Useful known results 

We now recall the Metzler matrices definition introduced by in mathematical economy. 
Definition 3 (Metzler Matrix) 

A matrix is Metzler iff all its off-diagonal elements are non-negative. 

We now focus on two results on Metzler matrices that are consequences of the original 
theorem of Perron- Frobenius on positive matrices The first comes from [llj : 

Theorem 2 

Consider a Metzler matrix A and a vector b £ int(R™ ); Then, A is Hurwitz iff: 

3x>0, Ax + 6 = 0. 

The second can be found in [20]: 
Theorem 3 

Consider a Metzler matrix A. Then: 



• The dominant eigenvalue of A (i.e. of largest real part) is real; its associated eigen- 
vector is positive. 

• consider a matrix B > A. Then the real part of the dominant eigenvalue of B is 
greater or equal to the dominant eigenvalue of A. 
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Let us recall the definition of the "cooperative systems", introduced by [2( (see also |21|) 
and that will be payed attention in the sequel. 

Definition 4 (Cooperative Systems) 

System £7J) is a cooperative system iff its Jacobian Df(x) is Metzler for all x € K™. 

One of the most interesting result on cooperative systems is to be found in |2lJ and reads: 



Theorem 4 

Consider a cooperative system 0) and two initial conditions xq and yo in M™ such that 
x o > Vo (resp. ^$>), Then: 

Vf > 0, x(t,x ) > y(t,y ) (resp. >). 



4.2 Preliminary result 

Proposition 2 (Positive Cooperative Systems) 

Consider a cooperative system Q). Then, the following are equivalent: 

• system yp is positive (resp. strongly positive) 

• /(0) > (resp. /(0) > 0) 



Proof 

Suppose /(0) > (resp. /(0) > 0). Let us denote the i-th row of Df(x) by Df t (x). Then 
we have for all X\ Xi =o > 0: 



Mx\ Xi=0 ) = mo) + 



Dfi(sx\ Xi = )ds 



X\ Xi =0- 



(7) 



Since Df(x) is Metzler for all x, the only possible negative term in the scalar product in {JJ 
is zero: indeed, Df iyi (the sole possible negative term of £>/») is multiplied by aij = 0. Then: 

Vx| X4=0 > 0, /i(x| X4=0 ) > M0) > (resp. > 0). 

Since this holds for all i £ [l--«], we conclude, using Proposition^! that the system is positive 
(resp. strongly positive). 



The reverse implication is obvious. □ 
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Remark 3 Notice that the class ofpositive cooperative systems is a generalization of the 
class of linear positive systems JTk . Q/. 

We now state a result on global asymptotic stability of an equilibrium point for a special 
class of strongly positive cooperative systems. 

Theorem 5 

Consider a strongly positive concave cooperative (SPCC) system fT|f. i.e. that verifies the 
following condition: 



Proof: 

We first show the exponential stability of a positive equilibrium x* (that is necessarily 
strongly positive) of a SPCC system Q. 

We obviously have: 



Since Df(x) is Metzler for all positive x, so does the bracketed matrix (denoted F(x*) in 
the sequel) in equation Since the considered system is a strongly positive cooperative 
system, we have from Proposition /(0) ^> 0. Applying Theorem we conclude that 
matrix F(x*) is Hurwitz. On the other hand, equation © implies: 



Since F(x*) is a Hurwitz Metzler matrix, we show using Theorem that its dominant 
eigenvalue is negative real. From the above inequation and Theorem we conclude that 
Df(x*)'s dominant eigenvalue is negative real too, proving that x* is exponentially stable. 

We now show the unicity of a strongly positive equilibrium x* . 

Suppose that the considered system possesses two strongly positive equilibria x\ and 
such that x* >x\. Then: 





Vse [0, 1], Df(sx*)>Df(x*). 



which yields: 




o = /«> 




(10) 
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From equation JH| and since x* > x 2 , we use similar arguments than in the previous part 
of the proof to show that: 



Df(sx\ + (1 — s)x 2 )ds 



<Df(x* 2 ). 



Since x 2 is a strongly positive equilibrium, Df{x^) is Hurwitz. Then, using Theorem^! the 
bracketed matrix is Hurwitz too. Then it possesses an inverse in A4 n (C) which, together 
with equation 11 Oil , implies that x\ = x 2 . 

Suppose now that the considered system possesses two strongly positive equilibria x* and 
x\, that are not linked by the relation ">", i.e. (x* — x%) is neither positive, nor negative. 
Let us define the parallelotopes B z C R™, such that: 

Vz g R" , B z = {x G R" , x < z) 

Using Theorem together with the considered system's cooperativity and strong positivity, 
it is easy to show that B x * and B x * are positively invariant sets. We now define the state 
£3, such that: 

x 3 , t = rmn(xi ti , x 2i ) 

Since B X3 = B x * (~l B x * , it is positively invariant by system £Q too. Then, using Brouwer 
fixed point Theorem (see e.g. 13]), there must exist a third equilibrium x\ belonging to 
B X3 . Note that, from x$ definition, we necessarily have: 

X3 < x\ and S3 < x 2 

This case has been treated previously and we conclude that what yields a 

contradiction. 

We conclude that if a SPCC system (QJ possesses a positive equilibrium, it is the sole 
positive equilibrium. 

We achieve the proof of Theorem0showing the global attractivity (on R" ) of the positive 
equilibrum x* of a SPCC system JQl. 

Consider the following real valued function: 

Vi G [l..n], Qi : R + — > R 

k ' * fi(kx*) 

Since /(0) 3> and f(x*) = 0, it is clear that <?i(0) > and <fa(l) = 0. Then, there exists a 
fco G (0, 1) such that: 

^(fc )= 5i (l)- 5j (0)<0. 
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From equation (jHJ, we get: 

0<k Q <k^0>^(k o )>^(k). 

showing that function gi strictly decreases from fco to +00. Since <?i(l) = 0, gi{k) < for all 
k > 1. This holds for all i <E [l..n], then: 

Vfc > 1, f(kx*) < 

Untill the end of the proof, we assume k > 1. Notice that, from the cooperativity and strong 
positivity properties, Bkx* is a positively invariant set. Now consider x, the time derivative 
of x. We have: 

x, = Df(x)x 

Since Df{x) is a Metzler matrix for all x, it is straightforward that if the initial state 
velocity x(t — 0) is positive (resp. negative) then it will remain positive (resp. negative) for 
all positive time. 

Remind that x(x — 0) 3> and x(x = kx*) <C 0. Then, the forward trajectory of system 
Q initiated at x — (resp. xq — kx*) is, in the ">" sense, an increasing (resp. decreasing) 
function of the time. Moreover, from B(kx*) positive invariance, it is lower (resp. upper) 
bounded by (resp. kx*), thus it converges, necessarily toward an equilibrium. 

x* is the sole equilibrium belonging to R™ , the two considered trajectories converge thus 
toward x*. From system cooperativity together with Theorem 0] we conclude that each 
trajectory initiated in Bkx* converges toward x* . Since this holds for any Bkx* (k > 1), we 
have shown that x* is globally attractive on K™ , which concludes the proof. □ 

Remark 4 Proof of theorem^ uses to a large extent the ideas proposed by Jst\] . However, 
his proof was dedicated to another class of systems: Kolmogoroff-type population dynamics 
models of cooperating species, which are based on different hypotheses. One can also have a 
look at 'lit] for some results related to Theorem^ 

Remark 5 It is clear that Theorem\^is quite similar to Hirsch's results /£/ on cooperative 
systems, rephrased by ft%] . and stating that for a bounded strongly monotone systemsfe.g. 
cooperative systems verifying xq > yo x(t,xg) 3> y(t,yo) , the trajectories, apart from a 
set of zero measure (thus at most n — 1 dimensional) of initial conditions, converge to the 
set of equilibrium points. Thus if such a system were bounded and the equilibrium unique 
then almost all trajectories converge to it. Notice also that Smith provides a related Theorem 
(see plj /. theorem 3.1). Here for SPCC systems, both uniqueness of the equilibrium (if it 
exists) and (real) global convergence are guaranteed, which is of particular interest for the 
application of our control procedure fjfi). 
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4.3 Sufficient conditions to verify HI 

On the basis of SPCC systems and Theorem0we provide sufficient conditions that are much 
simpler to check and that guarantees that Hypotheses HI are satisfied. 

Consider a system (TJJ under the following hypotheses: 
Hypotheses 2 (H2) 

1. Vx G int(K"), ip(x) > andVi, Ciip(x\ Xi=0 ) > 

2. /(0) > 

3. /(.) is such that system x = f(x) is cooperative 

4. Van, a* G R™ , x x < x 2 =► D/(a;i) > D/fe) 

5. 3/3 m G M+ smc/i V/3 > /3 TO , /3/(0) + c> 

6. V/3 > /3 m , 3a;* G int(R^) smc/i i/iai /3/(x*,) + c = 

One can show that H2 implies HI, thus that a system (2J under H2 is a candidate for 
Theorem ^ 

Indeed, H2-1 and Hl-1 are the same. H2-2 together with H2-3 implies through Propo- 
sition [2] that x — f(x) is a positive system i.e. Hl-2 holds. H2-3 together with H2-5 and 
Proposition implies that for all f3 > (3 m , the system: 

x = j3f(x) + c. (11) 

is a strongly positive cooperative system. This fact together with H2-4 implies that (1 1 H is 
a SPCC system. H2-6 ensures system ^] has, for all /3 > /3 m , an equilibrium x^ that is, 
through Theorem^ GAS on R^ so that Hl-3 holds. □ 

In the next section we provide numerical simulations, based on examples fulfilling hy- 
potheses H2, that illustrates Theorem interest and efficiency. 



5 Illustrative Examples 

We show on three examples that the possible open loop dynamics of systems J5J under H2 
might be very complex. Proofs of the complex dynamics are not detailed as they do not 
fall beyond the scope of this paper, but they can easily be performed with the help of the 
proposed references. 
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Example 1 (Bi-stability) This system models a simple bioreaction occuring in a contin- 
uous stirred tank reactor with substrate inhibitory effects (see e.g. JUj). This is a simple 
model of an anaerobic digester, an apparatus used for waste water treatment 0|. Both 
variables represent concentrations inside the reactor: x\ denotes the substrate (pollutant) 
consumed by the biomass (anaerobic microorganisms) xi to grow at a per capita rate fi(x2). 
u denotes the "dilution rate" (i.e. passing flow per volume unit) feeding the reactor with 
substrate at a concentration X\j n and withdrawing a blend of X\ and X2 from it. We get 
the following system that will be refered to as (Si) in the sequel. 



-1 
-1 

= ufi(x) + Ci1px(x). 

with: u(x{) — — ^ m 1 » 

Xi t i n , k, Hm, K m and Ki are positive constants. 

It is clear that Hypotheses H2-1, H2-2, H2-3 and H2-4 hold. Moreoever H2-5 holds with 

fim = k/xi,in- Now pick a j3 > f3 m , then (3fi(x) + c\ipi(x) — possesses a solution x*p such 
that: 

/ * 

p i » o. 



i 



V/3 > /3 m , x* p = 

\ 

so that H2-6 holds. Then we have through Theorem ^ 



Proposition 3 For all 7 > k/x\^ n the control law 

u(.)=7lMz) (12) 
globally stabilizes, on int(M"), system (Si) towards x* the (sole) solution of~ffi(x)+Ci = 0. 



Now let us consider the open loop behavior that might be produced by <Si. Suppose 
that Xi t i n > argmax and u £ (^(si !ln ), max(/i)) is constant. Then the considered system 
possesses three equilibria, two of which being stable on their respective basins of attraction 
separated by the stable manifold of the third equilibrium that is a saddle point. This 
behavior is depicted on figure^ see e.g. ^3] f° r a rigorous analysis. 

The two stable equilibria correspond for the strongly positive one to the "operating point" 
(the biomass survives inside the reactor and the pollutant concentration is reduced compared 
to its in the inflow) while the other do not (the biomass disappears and the waste water is 
no more treated). It is then necessary to ensure that no trajectory may be driven to this 
latter equilibrium. In this real life example, the quantity -01 (x) can be easily measured since 
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M Bi -stability behavior in open loop 




Figure 1: State space representation of various open loop forward trajectories of (S\) with o 
representing equilibria and u belonging to (/i(ii im ),max(/i)): existence of two locally stable 
equilibria. 

it corresponds to the outflow of biogaz that is produced by the anaerobic digester. Taking 
adavantage on this measurement, control procedure <PJ has been applied to a real life pilot 
scale anaerobic digester. It has shown its interest and efficiency for the global stabilization 
of such multistable devices. It thus avoids the crash of the system, see (lj| for more details. 

Example 2 (Attractive limit cycle) 
The following system is analogous to Goodwin metabolic models @. It will be refered to as 
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system (5 2 ) in the sequel. 



u 



fcl+Xl 



= uf 2 (x) + C 2 lp2{x) 




with I = 2.1, /zi = 2/2.1, /i 2 = 4 * (0.01 + 1/2.1), ki = 1/4.2, fc 2 = 0.01 + 1/2.1, a x = 0.01, 
a 2 = 1 - 2(0.01 + 1/2.1) and n = 80. It is easily proven that H2-1, H2-2, H2-3 and H2-4 hold 
true. Moreover H2-5 holds with f3 m = 0. Now pick a (3 > (3 m , then consider the equation 
Pf2{x) + C2ip2(x) — 0. We get that x 3 must be a solution of: 

MgQ^i +ai(fci/3i(l + q;g) + l)) 
X3 Q2 + M i + (ai+fc 2 )(fci^(l + x 3 ") + l) 

evaluating both sides of this equation at £3 = and +00, we show that there exist a solution 
x*p 3 which is positive for all /? > /3 m . x*p 2 and 1 are then computed and shown to be 
positive too, so that H2-6 is verified. 

Suppose that the considered system operates in open loop with the non-negative input 
u = 1. Then, as depicted on figure [2 (<S 2 ) possesses an attractive limit cycle around an 
unstable equilibrium point. The existence of a nontrivial periodic orbit can be proven with 
Theorem 1 from []}■ 

As illustrated on Figure we get from Theorem 



Proposition 4 For all 7 > the control law 

u(.)=yh(x) (13) 
globally stabilizes, on int(R"), system (<S 2 ) towards x* the (sole) solution o/7/ 2 (a;) + c 2 = 0. 

Example 3 (Chaos) The following system is a state space translation of the model of a 
chemical system so called "auto-catalator", see [ifj. It reads: 

7 "1 1 -1 k2 \ ( k ^-^W 

x =u — — \x+\ +1— \xixl, 

A0 1 \' -1) v k. )\ 

= uf 3 (x) + c 3 ip 3 (x). 
ki is positive, fc 2 E (0, 1) and k 3 > k± > . 
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Figure 2: State space representation of a forward trajectory of (52) with u = 1: existence 
of an attractive limit cycle. 



It is clear that H2-1, H2-2, H2-3 and H2-4 hold. Moreover H2-5 holds with (3 m = 

1/(^2(^3 — ^4))- Now pick a (3 > f3 m , then (3f 3 (x) + c 3 ip 3 (x) — possesses a solution x* } 
such that: 

/ [3k 2 k 3 + k 2 -l \ 



V/3>/3„ 



Xa = 




>0. 



since k 2 € (0, 1) and k 3 > k± > 0. Then H2-6 holds true. 
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Figure 3: Time varying representation of a forward trajectory of (S2) with u = 1 before 
t = 40 and under control law 11 -til from t = 40 to the end, with 7 = 2 > /3 m (= 0) computed 
s.t. asymptotically u(.) = 7^2(2;) equals its open loop value. 

(l^ | prove that Example El exhibits a chaotic behavior for the parameter values k\ = 
0.015, k 2 = 0.301, k 3 = 2.5 and u — 1 for all fc 4 . As a special case, we choose k± = 0.56. 
An example of a forward chaotic trajectory is shown on figure 21 

As illustrated on figure we get the following result from Theorem 

Proposition 5 For all 7 > the control law 

u(.)=jMx) (14) 
globally stabilizes, on int(M"), system (S3) towards x* the (sole) solution ofjfa(x)+C3 = 0. 
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Chaotic behavior in open loop 



X2 




06 1.5 X3 



Figure 4: State space representation of a forward chaotic trajectory of the system (53) with 
u = 1. 

6 Discussion 

In this paper we have considered a rather broad class of dynamical systems encompassing 
several possible dynamical behaviours. These models can especially represent biological 
systems exhibiting multistability, periodic solutions or even chaotic behaviour. In practice, 
the key point for stabilising these systems is that the adequate scalar function tp(x), whose 
exact analytical expression does not need to be known, is measured. In the example of 
the anaerobic digester that can have 3 steady states, the biogaz outflow rate ip( x ), ls easily 
on-line measured though its exact analytical expression depending on bacterial activity is 
rather uncertain. 

The proposed control law, although very simple, has proven its efficiency to globally 
asymptotically stabilize even chaotic systems which is known to be a rather tricky task 
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Stabilized chaotic forward orbit 
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Figure 5: Time varying representation of a forward trajectory of the chaotic system (£3) 
with u—1 before t = 20 and under control law itlljl from t — 20 to the end, with 7 w 1.73 > 
(3 m (~ 1.71) computed so that asymptotically u(.) = Jip3(x) equals its open loop value. 



when a part of the model (i.e. ip(x)) is not perfectly known. Moreover the control law 
naturally fulfils the input non negativity constraints which is classical e.g. for biological 
systems. 

The controller relevance was demonstrated in a real case for an anaerobic wastewater 
treatment plant (Mailleret et al., 2004). Despite its complexity (this ecosystem is based on 
more than 140 interacting bacterial species [3]) and the great uncertainties that characterize 
the biological models, the controller proved to efficiently stabilize the system. This approach 
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was also applied to a photobioreactor where microalgae consuming nitrate and light were 
growing [1J|. 

Another important issue when dealing with complex system is the difficulty to on-line 
measure the variable that must be regulated. The proposed control approach has thus the 
advantage of proposing a regulation scheme that does not require this measurement, which 
can be very difficult to carry out for biological systems. Moreover, and as it has been 
demonstrated with the real implementation, it is rather robust and stabilizes the system 
even if the model parameters are not perfectly known. The price to pay is however the 
possible lack of accuracy, since the solution of 7/(2;) + c = can differ from the expected 
set point in case of parametric uncertainty. In such a case the idea consists in adding an 
integrator to correct the bias. This adaptive approach was proposed in (Mailleret et al., 
2004) for the specific case of anaerobic wastewater treatment plant and the convergence of 
the controller was demonstrated. Of course it requires more measurement capability on the 
system, but we assume that the adaptation can be slower and thus needs a lower sampling 
frequency. The convergence proof of the adaptive controller in the general case is now the 
next challenge that we will tackle. 

Acknowledgements: The authors would like to thank V. Lemesle (COM - University 
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